arXiv:1508.04099v5 [quant-ph] 24 Apr 2017 


REVIEW ARTICLE 


Permanents, Bosons and Linear Optics 

Alexander Yu. Vlasov 

Federal Radiology Center (IRH) 

197101, Mira Street 8, Saint Petersburg, Russia 

E-mail: qubeat@mail.ru 

Abstract. Particular complexity of linear quantum optical networks is deserved 
recently certain attention due to possible implications for theory of quantum 
computation. Two relevant models of bosons are discussed in presented work. 
Symmetric product of Hilbert spaces produces rather abstract model. The second 
one is obtained by quantization of harmonic oscillator. In contrast to considered 
bosonic processes, so-called “fermionic linear optics” is effectively simulated on classical 
computer. The comparison of bosonic and fermionic case clarifies the controversy and 
the more elaborated oscillator model provides a deeper analogy. 
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1. Introduction 

Boson sampling |Tj is a formal model of nonuniversal quantum computations with 
noninteracting bosons suggested as a system that cannot be effectively simulated by 
classical computers. The model is much simpler than scalable universal quantum 
computer, but it might provide a principle evidence of so-called quantum supremacy 
important for theory of complexity. 

This paper is devoted to two issues. The first one is existence of a couple alternative 
quantum models relevant to description of the boson sampling. The second question 
is comparison with “fermionic” quantum circuits effectively simulated by classical 
computer. 

In section[2]the “ abstract” bosons model is briefly discussed and after short reminder 
about linear optics in section [3] more elaborated oscillator model is discussed in section|4j 
The comparison of “bosonic” and “fermionic” case is provided in section [5j 

2. Abstract bosons model 

An abstract quantum system with d states is described by d-dimensional Hilbert 
space %. An arbitrary state can be expressed using Dirac notation [2] with a basis 


Permanents, Bosons and Linear Optics 


2 


d d 

\^) = J2 c k\k), j>i a = i- ( x ) 

k =1 /c=l 

Group U(rf) of unitary matrices describes an evolution of the state. In such 
description a global phase does not matter [2] and the group SU (d) of unitary matrices 
with unit determinant may be used instead of U(rf) 

d 

I VO = u\4>), c'j = Yu jk c k , u G SU(d). (2) 

k =1 

Let us denote S n [PL) linear space of symmetric n-tensors, dim S n (PL) = C^ +n_1 . 
Symmetric product of tensors also can be used with similar purpose. The method is 
related to linear space of polynomials of degree n with d variables |3]. For n — 1 
such a space is simply PL. The S n (PL ) can be considered as a state of system with n 
indistinguishable bosons [3]. An element of S n (PL ) represented as symmetric product 
of n vectors from PL is called decomposable. It may be treated as a system with n 
noninteracting bosons. 

Two methods of indexing are used for basis of S n {PL ) pE|. The first one is sequence 
j with n indexes 

3 = 1 < ji < J2 < < j n < d. (3) 


Such notation is obtained from initial basis of tensor product after symmetrization of 
indexes. It may be rewritten due to polynomial representation [1] 


Xj 1 X j 2 


■ X jn — 


a . • • x/, rf + r\ H-1- r J d = n, (4) 

where r J k > 0 is number of indexes k in the sequence j. The numbers r J k produce an 
alternative representation as a vector with d elements. The notation is convenient 
for expression of standard normalization for basic vectors in space of symmetric tensors 
or polynomials 


T s _ 
x j 


./■ 


x„ 


yrj 


Tj — r\ \... r 


( 5 ) 


-i 


and the element of S n {PL) may be considered as a complex vector with C^ +n 
components in a basis | Xj). Action of unitary group on each PL for such normalization 
corresponds to an unitary representation on the S n (PL ) with respect to usual Hermitian 
scalar product and the unitarity [4j is natural due to 


^ ^ x k x k 

k =1 


52 r k= n 


n] n x k x k 


r i! 


■ r d ! 


= n 


( 6 ) 


jes 


For already mentioned earlier decomposable elements of S n (PL ) relations between 
the scalar products on S n (PL ) and PL may be expressed using permanents [5j. Let us 
recall, that the permanent of n x n matrix A is 

n 

per (A) = Aj,a(j), 

a j =1 


( 7 ) 
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where cr denotes all possible permutations of indexes. The only difference with 
determinant is lack of minus signs for some terms. However, unlike the determinant 
the permanent is an example of computationally hard problem [6j. 

Transformation of product (j4j) for change of variables x'■ = Ylt=i Mjk x k can be 
expressed with some polynomials V{M) of degree n with elements of the matrix M 

( 8 ) 


x 'h x 'h 


■ x „• = 




x ki x k-2 • • • x k n 


fci, 


A simple case suitable for further examples is n = d with consequent indexes k — j — 
(1, 2,..., n), when the polynomial simply equal to permanent 

= ( 9 ) 


More general case should be discussed elsewhere and analogues of (J9]) also could be 
expressed using permanents of matrices composed from rows and columns of M. 

An application of such an abstract model for discussion about permanent 
and determinant complexity may be found in [7]J together with suggestion about 
ineffectiveness of finding permanents using quantum processes with bosons because of 
big variance of measurement outcomes. 


3. Linear optics 

The abstract bosons model described above says a little about physical realization. 
Alternative way of producing expressions related with permanent may be obtained using 
quantum model of linear optical networks [8]. It may have relevance with theory of 
quantum computing [9] and computational complexity of linear optics was discussed 
further coining the term boson sampling [T] stimulating the series of experiments 

mammalia. 

An evident distinction in formulation of such a model in comparison with the 
abstract bosons model from section [2] is definition of some basic transformations using 
creation and annihilation operators dj, d), j = l,...,d (HUH]. 

The approach is quite natural, because with such a notation linear optical network 
with conserved total photon number corresponds to transformation [9] 

d 

a, = y Ujkdk, U G U(d), j,k — l,...,d. (10) 

fc=i 

Due to some analogy between (fTOT) and ([2]) they might cause similarity in formal 
calculations, but meanings of the expressions are quite different. The (1TO describes 
transformation of operators in Heisenberg representation, but (J2|) is applied to states. 

Transformation of some operator A due to (EJ) in the abstract bosons model from 
section [2] in Heisenberg representation ci would be 

A' = uAu ] . (11) 

Here u produces the same result as e^u reaffirming sufficiency of using u G SU (d). In 
(HU such compensation of phases is dispensable and so the whole unitary group U(d) 
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may be implemented, but it is rather a hint on more essential difference between the 
models discussed further. 

In fact, a model with symmetrization of states very similar with abstract bosons 
model from section [2] with infinite-dimensional space PL sometimes used in quantum 
optics as well na. It can be asked in turn, how to rewrite (HU in a way similar 
with CD- For such a purpose in next section transformation with conservation of 
total photon number is considered as particular case of most general linear Bogoliubov 
transformations pun without such requirement. 


4. Oscillator model 

4-1. Schrodinger description 

The model of quantum harmonic oscillator is recollected below with a brief excursus 
into theory of symplectic and metaplectic groups necessary for applications to linear 
optics pH El HE]. Let us consider infinite-dimensional Hilbert space PL and operators 
q, p of coordinate and momentum. In Schrodinger description the PL is associated with 
space of wave functions G L 2 (M) (square integrable) and the operators q, p are 
defined as 

q: -0 (q) —> qfi>(q), P' if{q) — > ~ idip(q)/dq (12) 

with canonical commutation relation (CCR) 

[p, q]=pq-qp = -i, (13) 

where system of units with h — 1 is used for simplicity. The generalization on set of 
operators qk, fik, k — 1,..., d with CCR 

\pj,Qk] = j,k = l,...,d (14) 

is straightforward using space L 2 (R d ) of wave functions with d variables. 


4-2. Symplectic group 

A real 2d x 2d matrix A preserving bilinear form 

d 

(*£j ?/)a ^ ^ fi'kljd.+k Xd+kVk) 

k =1 

is called symplectic. It also has property [19] 

A T JA = J, 

where A T is transposed matrix and J is 2d x 2d matrix 



(15) 


(16) 


( 17 ) 
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with O t / and l ( i are zero and unit d x d matrices. A composition of such matrices also 
satisfies (fT6|) and, so, the symplectic group Sp(2d, R) [T6j (TTj is defineclj] in such a way. 

Due to an analogy with orthogonal matrices R G SO(d) satisfying R T R = 1 and 
preserving Euclidean norm (x,i/)e = Yh x kVk a matrix with property (TT6l) and arbitrary 
J sometimes is called J-orthogonal. Such a definition includes both orthogonal (for unit 
J) and symplectic matrices [19]. 

Matrices preserving both symplectic (jT5]l and Euclidean forms belong to orthogonal 
symplectic group defined as intersection Sp(2d, R) 0 SO(2d). The group is isomorphic 
with unitary group U(d) [19] . Indeed, if to consider complex variables 

Zk %k T Dfc+dj Zk — •Kk Dfc+d; (1®) 

the real and imaginary part of Hermitian complex scalar product correspond to 
Euclidean and symplectic forms, respectively (z\z') = (x,x')e + fix,x')\. 

To show relation of symplectic group to CCR let us write fi ,, fik as formal vector of 
operators with 2d elements [T6j 

w = ( w k ) = (<7i, ■ ■ ■, QchPu ■ ■ ■ ,Pd)- (19) 

Equation (El) can be rewritten in such a case as 

[wj,w k \ =-iJ jk , j, k = 1,... , 2d, (20) 

where J 3 k are elements of matrix J (TTTTl . Due to such property 2d operators w' 

2d 

Wj = ^2 S ikWk, j,k = 1,..., 2d, (21) 

k=1 

also satisfy (|20T) if matrix S G Sp(2d, R). 

4-3. Metaplectic group 

Both the sets of operators up and wl related by f[2Tj) satisfy some form of CCR (I20T) 
and in agreement with general results about uniqueness of CCR they should be unitary 
equivalent , i.e., for any matrix S in (l2Tf some unitary operator Us should provide 
transformation mmm 

w'j = UsWjUg 1 , j — l,... ,2d. (22) 

Due to (1221) Us and e^Us correspond to the same matrix S, but such a phase freedom 
may be withdrawn and the only inevitable ambiguity is a sign ±Us- The group producing 
such a 2-1 homomorphism on Sp(2d, M) is known as metaplectic, Mp(2d, R) [HJCETls 18] - 
The unitary representation of Mp(2d, R) used in (l22li is not a finite-dimensional matrix 
group, but can be expressed by exponents with appropriate linear combinations of WjWk- 
The symplectic group is relevant also to classical optics, but metaplectic group is 
essential in quantum case [Ml El HE]. The group Mp(2d, R) describes most general linear 
optical networks with d modes and a subgroup of transformations with conservation of 
total photon number is discussed below. 

t Slightly different parametrization Sp(d, R) is used in [181119] . 
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The transformations respecting also Euclidean norm, i.e., sum of Hamiltonians of 
harmonic oscillators 

..2 d d ^ 

H=2 = H t = -(pl + q 2 k) (23) 

k= 1 k =1 


correspond to already mentioned orthogonal symplectic group isomorphic with unitary 
group U(d). Complex coordinates (TT8l) now associated with annihilation and creation 
(“ladder”) operators 


1 


a k = ^={q k + i p k ), a\ = ~^(q, 


ipfcj 


(24) 


v'T n " V2 

Equations (ITU [2~T1) are widely accepted Usmans], consistent with rather standard 
map to U(d) [T9] and used further in this work for simplicity instead of more general 
versions also relevant to quantum optics 
1 


H=-(p 2 + u 2 f), 

a = , (c oq + ip), a 1 = (uq — ip). 


(25) 


\/2cJ 


\f%j 


The oscillator model relates group U(rf) treated as a subgroup of Sp(2d, R) with 
conservation of “total photon number” defined by operator 

d 

* = £ 4 "i = ( 26 ) 

3 = 1 

Instead of f[2T]) already mentioned earlier should be used (TTOl) and expression for a) 
is obtained using Hermitian conjugation 
d d 

4? = E u iA, U€\j(d). (27) 

k =1 k= 1 

Equation (OTj) rewritten in such a way without requirement about photon number 
conservation would include terms a', ci in both sums (1271) reproducing most general 
Bogoliubov transformations. 

In such an approach subgroup of Mp(2d, M) corresponding to U(d) sometimes is 
denoted as MU (d) [2T]J and can be expressed by exponents with linear combinations of 
ajafc. An analogue of (l22j) is 

a) = LlijdjLtj) 1 , j — 1,..., d, W 6 MU(d). (28) 

with conjugated expression for aj. 

The MU(cf) is double cover of U(d) with sign ambiguity inherited from relation 
between Mp(2d, M) and Sp(2d, M) [18], 21]. The additional subtle problems may appear, 
because using some formal manipulations with lAu an one-to-one map with U(d) could 
be obtained, but it cannot be extended to the whole Sp(2d, M) [18]. 
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The difficulties with unitarity after application of discussed trick to get rid of 
double cover would produce rather unrealistic model of photons with prohibition to 
use squeezing transformations. Even if requirement of particle number conservation is 
justified for model with a massive bosons the approach discussed in this section has 
other important differences from the abstract bosons model introduced in section [2} 

A noticeable formal distinction is an additional phase parameter, because action of 
U (d) is not reduced to SU(d). The nontrivial structure of Ur is manifested here, because 
(1281) is not sensitive to formal phase multiplier C^IAr. So action of phase multiplier 
on matrix U in (1271) is implemented in alternative way and “embedded” directly into 
structure of U via additional term proportional to total photon number operator N (1261) 
in the exponent with quadratic expressions for elements of group mentioned earlier. 

Yet another important property of used model with double cover and sign ambiguity 
due to (j22]h (T28l) is similarity with relation between orthogonal and Spin groups [18, 21]. 

5. Comparison with formal fermionic model 

The theory of Spin groups is relevant with question about difference in complexity 
between bosons and fermionic model often associated with matchgate circuits effectively 
simulated on classical computer. The matchgate model was introduced in [22] with 
reformulation using formal fermionic model in [23} [ 21 ]. Similar fermionic model and 
equivalent approach with Spin group and Clifford algebras was also applied earlier to 
quantum computing problems in [25} [26} (27} [28]. 

The theme was further developed in works devoted to effective classical simulation 
of such a class of quantum circuits and they are often described using fermionic 
annihilation and creation operators [29} [30, [31]. The notation fij is used here for 
such operators for distinction from bosonic case. The analogues of uij (TTU]1 here are 2d 
generators of Clifford algebra 

hj -i = a,- + fij, c 2j = -i(a,- - ht). (29) 

Transformation properties operators (1291) are similar with (I2T[) . (!22[) and may be written 
as 

2d 

S R CjSf l = ^2R jk c k , (30) 

k =1 

but Rjk now are elements of orthogonal matrix and Sr G Spin(2d) is corresponding to 
matchgates (or “fermionic”) quantum circuit. The Sr has matrix representation, but 
number of components is exponentially bigger than in matrix R , because it corresponds 
to quantum network with d qubits. 

A restricted case of transformation conserving number of fermions discussed in 
[23] provides even closer analogy with boson case. If the evolution is expressed as 
exponent of Hamiltonian with linear combination of fijfi fc then the fermionic operators 
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are transformed [23] by direct analogue of (12711 

d 

Suq} 3 Su — ^ Ujk^ki ( 21 ) 

k =1 

with Su is from subgroup of Spin(2<i) corresponding to restriction of SO (2d) on 
“realification” of SU(d), i.e., orthogonal symplectic group mentioned earlier. 


Two cases should be taken into account for comparison of complexity fermionic and 
bosonic models. The single-mode measurement in terminology of [29] is the first one 
and mainly used in [30( j31], 32]. In the concise form the setup [31] for such experiment 
is any computational basis state |xi... Xd) as the input and a final measurement of an 
arbitrary qubit in the computational basis as the output. For arbitrary state | if) such a 
measurement may be described by probability p\ ^ to obtain 1 as a result of measurement 
for a qubit with index k 

Pi ] = (V’lOfcOfclV’) (32) 


and an analogue of such expression in bosonic case 

= (T|7V fc |T) = (N k ) (33) 

is corresponding to expectation value for number of particles in a mode k. After 
application of linear optical network new expectation values (N' k ) are 

i*')=mi« , >. <K) = wyuviw 

It can be rewritten using (1271) and (1281) with matrix U corresponding to operator U 


(34) 


iK) = 0 (E M) (E u k & 

3 = 1 j =1 

d 

= E 




(35) 


3 , 3 '= 1 


If an effective way to calculate (TlatdjvlT) exists for given input state IT) then the 
(N' k ) also can be calculated in poly time using (135|) and the methods are similar with 
matchgate model mm- 

The Fock states [15j can be considered as an alternative for computational basis for 
input state for a case with n bosons and d modes 


|Tj) = 




1 < jl ■ ■ ■ < jn < d 


(36) 


with |0) is vacuum state. The coefficients a/T^ can be derived from expressions for 
ladder operators |T5] and coincide with normalization (]5]) for vector = {r{,... ,r 3 d ) 
defined by (J4J) for given sequence j = (ji, J 2 , • • •, j n ) introduced in (J3J). So, finally 

(zUP iq\ 
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For such states (\& J -|a£afc/|\E r j) may be effectively calculated. 


The other case is multi-mode measurement [23j, [29]. The analogues of such model 
in bosonic case would use instead of (1331) polynomials with d' k , a k of higher degree. 
So, different computational complexity of determinant and permanent may be really 
essential in some examples. 

It may be written U |0) = 9^ |0) with some phase multiplier = 1, because U 
is expressed as an exponent with linear combinations of djdk and for vacuum state 
dj |0) = 0. Therefore, 


U\Vj) 


-^Ua)d] 2 ---d)u 1 W|0) 


(38) 


After further application to each multiplier 1Aa}-U~ x and expansion as sums using (1271) 
expressions with permanents may be obtained. 

Anyway, ()33l) and analogues with small number of ladder operators may hide 
complexity arising from application of linear optical network to Fock states. An example 
is (134|) or (l35j) that may be used for effective classical computation of single-mode 
measurement outcomes (expected average number of photons ( N ' k ) for each mode k ) 
after application of linear optical networks to Fock states. 

The permanent complexity is more essential in expressions for “transition” 
amplitudes (and probabilities) between two Fock states 


= (d'fc|W|'F ? -), p kj = \a kj \ 2 . (39) 

Basic example with consequent indexes in k and j corresponds to ([9]). For 
calculation of probabilities pkj a phase |6^| = 1 from (|38l) should be omitted. An 
alternative consideration with commutative polynomials resembling an abstract boson 
model discussed earlier in section M may be found in [U 133] . 

Analogue models with preserving number of fermions was studied in [23]. The same 
d- dimensional unitary group 651 again can handle the evolution, but determinants 
are used instead of permanents and expressions for fermionic amplitudes may be 
efficiently evaluated [23]. The analogue of (1391) in fermionic case directly coincides 
with determinant. So, for fermions multi-mode measurement amplitudes also can be 
efficiently evaluated due to absence of problems with permanent calculation discussed 
earlier for bosons. 


6. Conclusion 

Some topics relevant to consideration of complexity of simulation of quantum processes 
with boson were discussed in this paper. A distinction between an abstract bosons model 
and more elaborated approach with quantum harmonic oscillator was emphasized. The 
treatment of real photonic system may require consideration of even subtler problems 
that should be discussed elsewhere. 
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Let us recollect similar equations for transformations of some sets of states | ipk) and 
operators Cjj relevant to many examples discussed above using notation 

Wk) = S\i>k), O'j = SQjjS~ l = 22 SjkOk- (40) 

k 

For oscillator model of general linear quantum optical network operators u 3 = vbj are 
defined in (HU and the Sjk is symplectic matrix (j2T|) with S is unitary operator Us (l22l) 
from metaplectic group. For an important case with conserving number of bosons Cjj 
correspond to annihilation and creation operators, S was denoted as Uu in ([28]) . and 
Sjk is unitary matrix Ujk in (j27l) . 

In section [5] about fermionic model Cjj = c 3 (129|) . S = Sr G Spin(2d) and Sjk 
corresponds to orthogonal matrix R (1301) and conserving number of fermions is taken 
into account in (13T|) with unitary matrix. 

Therefore, for conserving number of boson and fermions evolution can be described 
by the same unitary group with “reduced” dimension d instead of 2 d for most general 
quantum circuit and significant difference in complexity of amplitudes evaluation for 
multi-mode measurements looks especially challenging. 

The fermionic model may be effectively simulated by classical computer producing 
some controversy with bosonic case. On the other hand average photon numbers in 
output of each mode can be simply calculated likewise with fermionic case. So, classical 
computer could effectively simulate output for each mode, but without proper quantum 
multi-photon correlations. Recent achievements in experiments with photons [33, 31] 
permit to check such subtleties. 

For abstract bosons model introduced in section [2] unitary operator S acting on 
states could be compared with u G SU(<i) in (J2J) , but without analogues of operators 
such as Cjj in (1401) . 

A Fock state (136|) may be considered as an analogue of ([5]) in abstract bosons model , 
but confusing them may produce certain problems. Let us consider simplest case with 
action of Ujj on a basis of PI 

|fc) = a*|0). u v \k) = Z44l0> = zya'wpz^iB) 

= E UiSfiaW = e uT, u h\ k > = 6 uU’\k), (41) 

3 3 

where \6^\ = 1 was already introduced earlier, |0) denotes vacuum state, and | k), 
k = 1,..., d correspond to basis of PL. In such a case in abstract bosons model difference 
between U\j and U* lacks of proper treatment and simply missed sometimes. On the 
over hand (127|) and (1281) 

d 

a'j = 22 Ujkdk = UudjU-u 1 . (42) 

k= 1 

would not make a sense, if Uu and U* are not distinguished. Indeed, both signs ±U 
lead to the same matrix U in (l42j) and define of Uu as an operator from double cover of 
unitary group. 
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Together with the necessity in consistent mathematical expressions there are also 
physical reasons to make distinction between two models discussed above. The oscillator 
model based on quantization of classical linear optics and the abstract bosons model is 
close related with quantum approach to discrete models such as group of permutations 
[ 4 ]. In such a way the subtle relations between two models may be compared with a 
wave-particle duality. 
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